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ABSTRACT 


The  dynamic  stability  of  a  parachute -load  system 
has  been  analytically  investigated  for  a  point-mass  load 
and  a  statically  stable  parachute.  A  typical  system  con¬ 
sisting  of  a  relatively  large  suspended  load  mass  and  small 
ribless  guide  surface  parachute  has  been  numerically  cal¬ 
culated.  Utilizing  the  apparent  mass  and  apparent  moment 
of  inertia,  as  well  as  the  aerodynamic  coefficients  of  the 
parachute  canopy,  the  equations  of  motion  for  the  system 
have  been  solved.  The  influence  of  several  design  parameters 
upon  the  dynamic  stability  characteristics  of  the  system 
has  been  discussed. 
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angular  velocity  of  the  parachute  axis 


I.  INTRODUCTION 


When  considering  the  complete  pattern  of  motion 
of  a  parachute -load  system,  the  dynamic  stability  of  the 
system  is  a  prime  consideration. 

For  such  a  consideration,  one  thinks  of  a  freely 
moving  system  as  having  six  degrees  of  freedom,  consisting 
of  linear  and  angular  velocities  about  three  mutually  per¬ 
pendicular  axes.  The  angular  velocities  result  in  an 
oscillating  motion,  with  the  mode  of  oscillation  generally 
referred  to  as  the  stability  characteristic  of  the  system. 

It  has  been  noted  (Ref  l)  that  such  characteristic 
motions  depend  on  the  aerodynamic  coefficients,  their  deriva¬ 
tives,  and  on  the  mass  and  moments  of  inertia  of  the  mechan¬ 
ical  system  as  well  as  on  the  apparent  mass  and  apparent 
moment  of  inertia. 

Combining  an  aerodynamically  neutral  load  with  a 
static  aerodynamically  stable  parachute,  a  closed  form 
solution  is  possible  provided  that  the  initial  oscillations 
are  small  and  that  the  parachute  has  a  constant  stability 

Qc 

derivative  ( - — )  in  the  range  of  oscillations.  This 

vacc'3 

process  is  valuable  for  a  certain  group  of  applications 
and  has  the  advantage  of  being  acceptable  for  engineering 
calculations . 

For  a  more  general  type  of  problem,  where  the 
oscillations  may  be  large  and  the  stability  derivative 
is  not  constant,  solutions  can  be  obtained  merely  for  specific 
cases  and  in  a  numerical  way  by  means  of  analog  or  digital 
computers.  Several  of  these  more  general  as  well  as  specific 
cases  have  been  investigated  and  solved  in  a  series  of  pub¬ 
lications  by  R.  Ludwig  and  W,  Heins  (Refs  2,  3  and  4). 

These  studies  discuss  also  the  influence  of  several  design 
parameters  such  as  suspension  line  length,  effective  porosity 
of  the  canopy,  apparent  mass  and  weight  of  the  load  for  the 
particular  application. 


Manuscript  released  by  authors  February  1964  for  publication 
as  a  FDL  Technical  Documentary  Report. 
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II.  EQUATIONS  OP  MOTION 


For  the  analysis  of  the  dynamics  of  a  parachute - 
load  system,  one  assumes  the  case  of  two  falling  bodies 
connected  by  rigid  lines.  Figure  1  presents  the  load- 
parachute  system,  the  body  fixed  coordinate  axes,  x-y-z, 
as  well  as  the  system  of  external  forces.  The  x-axis  coin¬ 
cides  with  the  parachute  axis  of  symmetry,  the  y-axis  is 
perpendicular  to  the  x-axis  and  in  the  plane  of  motion  of 
the  system  with  the  z-axis  defined  by  a  right  hand  coordinate 
system.  The  origin,  (0,0,0),  is  at  the  center  of  mass  of 
the  entire  system. 

In  deriving  the  equations  of  motion  of  the  system, 
the  following  assumptions  will  be  used: 

1.  The  entire  system  constitutes  a  rigid  body. 

2.  The  load  is  considered  as  a  point  mass  and 
does  not  possess  aerodynamic  characteristics 
or  a  moment  of  inertia. 

3.  The  mass  and  the  aerodynamic  forces  of  the 
suspension  lines  can  be  neglected. 

4.  The  effect  of  the  apparent  mass  acts  at  the 
center  of  volume  of  the  canopy. 

5.  Ths  motion  is  restricted  to  the  x-y  plane. 

The  velocity  of  the  center  of  mass  of  the  entire 
system  can  be  presented  as : 


v  =  vxi  +  vy5  (1) 

where:  vx  =  velocity  of  the  center  of  mass  in  the 
direction  of  the  parachute  axis 
Vy  *s  velocity  of  the  center  of  mass  in  the 

direction  perpendicular  to  the  parachute  axis. 

The  corresponding  momentum,  M,  of  the  parachute -load  system, 
excluding  the  effects  of  the  apparent  and  enclosed  mass,  is: 

H  =  (mp  +  mj  )  (vx£  +  vy5)  (2) 

where:  =  mass  of  parachute  material 

mjj  =  mass  of  the  suspended  load. 
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The  equations  of  motion  of  the  system.,  lateral 
and  rotational,  can  be  written  in  accordance  with  Newton1 
Law  as : 


a) 


_  d^1)  M 

^-F  “  - 


b) 


H 


s 


(3) 


in  which  the  'following  notation  is  used: 

P  =  external  force 
R  =  radius  vector 

angular  momentum. 


H  = 

The  derivative 


dt 


LI 


represents  differentiation 


with  respect  to  a  space  tfixed  (inertial)  reference  frame 
(x 1 ,  y',  z')  (Pig  2).  A  derivative  with  respect  to  a  par¬ 
ticular  reference  frame  can  be  expressed  with  respect  to 
another  frame  which  rotates  relative  to  the  first  one  by 
means  of  the  relationship  (Ref  5,  pp.  53-55): 


where : 


+  lo  x  A 


A 

cTE 

d<S>A 

cTE 


rate  of  change  of  A  with  respect  to 
reference  frame  1 

rate  of  change  of  A  with  respect  to 
reference  frame  2 


(4) 


co  =  angular  velocity  of  reference  frame  2 
with  respect  to  reference  frame  1. 


In  this  case,  A  =  Iff,  reference  frame  1  is  the 
space  (x',  y’,  z')  frame,  and  reference  frame  2  is  the  body 
fixed  reference  frame  (x,  y,  z). 


Therefore,  one  may  write: 


+  co  x  M 


(5) 


where  uj  =  oj  k  =  CaJ  k!  =  angular  velocity  of  the  para¬ 
chute  axis  with  respect  to  the  z-axis. 


Performing  the  operations  indicated  in  this  equa¬ 
tion,  one  obtains: 
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M 

cTE 

Pig  1: 


(mp  +  irijj )  [(vx  -  vycJ)i  +  (vy  +  vxo;)3]  .  (6) 

The  aerodynamic  forces  arej  in  accordance  to 

*a  =  -Ti  -  Nj  (7) 


and  the  gravity  forces  are  given  by: 

Pg  =  (Wj i  +  Wp)  cos  0  i  -  (W£  +  Wp)  sin  ©  5.  (8) 

Note:  the  gravity  force  on  the  enclosed  mass  is 

balanced  by  the  bouyancy  force  on  the  canopy. 

The  effects  of  the  apparent  and  the  enclosed 
mass  are  given  by: 


Fam  -  -  ”ayayoy?  (9) 

where:  m„  =  apparent  mass*  in  the  x-direction 

ax 

ma  =  apparent  mass*  in  the  y-direction 

y 

ax  =  acceleration  of  the  center  of  volume  of 
cy-  the  canopy  in  the  x-direction 


cv. 


=  acceleration  of  the  center  of  volume  of 
the  canopy  in  the  y-direction. 


The  acceleration  of  the  center  of  volume  of  the 
canopy  is  given  by: 


c.v. 


“ar 


—  —  3^;  1  "t  9-y  J 

C.V.  ^  C.V. 


(10) 


where.:  vcv 


velocity  of  the  center  of  volume  of  the 
canopy  with  respect  to  the  space  fixed  frame 


vxi  +  vy3  +  x(-L2i) 

vxi  +  (vy  -  co  L2)?. 


(11) 


Therefore j  utilizing  relation  4  one  obtains  the  total 
acceleration: 

*See  definitions  of  symbols . 
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ac.v.  =  [vx  -<*)  (vy  -  CJ  Lg)  ]  i  +  [vy  -  ^  L2  +  w  vx]  3  .  (12) 

Consequently,  the  acceleration  In  the  x-  and  y-directions  Is: 

=  vx  -  (vy  -  uj  L2) 


cy 


’Vc.v.  =  ^  l2  +  w  vx 


(13) 


Utilizing  these  relations  in  Eqn  9»  one  finds  the  apparent 
mass  force  as: 

Pam  =  -max[vx  -wvy  +  U)2  L2]  i  -  may[vy  -  tb  L2  +^vx]  J.  (14) 

If  one  now  substitutes  all  of  these  forces.,  and 
relation  6  into  Newton's  Law  (Eqn  3a),  two  scalar  equations 
representing  the  motion  in  the  x-  and  y-  directions  are 
obtained. 

(nip  +  m&x+  m£)(vx-vy^)+  maxL2  +  T  -  (W^  +Wp)cosG  =  0  (15) 

(nip  +  ma  4-m^ )  (vy+vx )  -  mayL2k>  +  N  +  (W^  +Wp)sin  0  =  0.  (l6) 

y 

Newton's  second  law,  which  governs  the  rotational 
motion  (Eqn  3b),  can  be  expressed  as  follows: 


d  ^  —  =  [Rxf'  (17) 

with: 

H  =  (lom  +  Ia)uk  (18) 


where :  H 


R 


F 


angular  momentum  of  the  system  with  respect 
to  a  nonrotating  reference  frame  with  its 
origin  at  the  center  of  mass  of  the  system 

position  vector  from  the  center  of  mass  of 
the  system  to  the  point  of  application  of  an 
external  force  F 

external  force 


Ia  =  apparent  moment  of  inertia  of  the  parachute 

canopy  and  inertia  effects  of  the  enclosed  air 
about  the  center  of  mass  of  the  system 
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^cm  =  rnoment  of  inertia  of  the  load  and  parachute 
material  about  the  center  of  mass  of  the 
system 

The  external  moments  (  L*  x  7  )  are  given  by: 

x  F  =  (-L1i)x(-Nj)  +  )x  (-Wjg  sin  Q  3)+( )x ( -WpSin  9  j) 
or 

x  f  =  k  [+  L^N  -  W5L3s  in©  +  WpL^sin©] 

and  since  the  last  two  terms  [-  W^L^sin  e  +  WpL^sine]  cancel 

each  other  because  of  the  definition  of  the  center  of  mass, 
one  obtains : 

[r  x  P  =  kL-jN  (19) 

The  moments  due  to  the  apparent  mass  effect  which, 
in  this  discussion,  encompasses  the  effect  of  the  enclosed 
air  mass  have  not  been  identified  separately  since  their 
contribution  to  the  moment  equation  is  included  in  the 
experimentally  determined  value  of  I  . 

a 

Introducing  relations  18  and  19  into  Eqn  17  yields: 

(l)  f  ]  a 

V-  [dcm  +  =  +  L1N*  (2°) 

and  after  rearranging: 

km  +  ij  «  -  L1N  "  0  '  <21> 

Three  differential  equations  for  the  motion  of  the 
parachute -load  system  have  now  been  obtained.  They  are  re¬ 
produced  here  for  future  reference: 

momentum  equation  for  the  x -direction 

(mp+max+m^  )  (^x-vy^)+  m^LgOJ2  +  T  -  (W,  +Wp)cos  ©  =  o  (22) 
momentum  equation  in  the  y-direction 

(iT‘p+rnay+m  i )  (vy+vx  to )  -  mayL2^+  N  +  (W£  +Wp)sin  ©  =0  (23) 

angular  momentum  equation 

dcm  +  Ia)  ci  -  LxN  =  0  .  (24) 
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III.  LINEARIZED  EQUATIONS 


Equations  22  through  24  constitute  the  general 
equations  of  motion  of  a  parachute -load  system  and  apply 
equally  well  to  large  or  small  oscillations.  It  should 
he  noted,  however,  that  they  are  nonlinear,  i.e.,  con¬ 
taining  terms  of  the  form  vy  uj  .  Thus,  it  becomes  ex¬ 
ceedingly  difficult,  if  not  impossible,  to  find  a  closed 
form  solution  of  these  equations  without  reasonable  restric¬ 
tive  assumptions. 


The  simplest,  but  perhaps  also  the  most  important 
case  of  the  dynamic  stability  problem,  is  given  by  a 
parachute -load  system  in  which  a  parachute,  which  is  stati¬ 
cally  stable  about  zero  angle  of  attack,  oscillates  over 
a  relatively  small  angle  range  whereby  its  stability  de- 


assumes  a  constant  value. 


Furthermore , 


on  ^R’ime  that  the  system  descends  approximately 

vertically  v.^h  its  equilibrium  velocity  where  the  aerodynamic 
drag  equals  the  weight.  These  assumptions  may  be  expressed 
as : 


sin 

9 

=  e 

.  cos 

e 

1 

sin 

a 

=  a 

cos 

a 

1 

sin 

A 

cos 

A  = 

1 

r  ^ 

v  = 

V 

o  +  v 

L v 

«  voJ 

where  vQ  =  equilibrium  velocity 

'v  =  deviation  of  v  from  the  equilibrium  velocity. 


From  the  geometry  of  Fig  1: 


vx  = 

v  cos  a  =  (v0  +  ’ 

v)cos  a 

ii 

> 

-v  s  in  oc  =  -  (v  + 

v)sinO( 

vx  = 

-(v0 

+  ~ )  s  in  oc  oc 

+  v  cos  a 

vy  = 

-K) 

+  v)cos  oca 

-  v  sin  a 

(26) 


Using  the  small  angle  assumptions,  and  neglecting  second 
order  terms,  i.e.,  terms  of  the  form  a  a  ,  one  obtains: 
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(26a) 


'y 


vQ  +  v 

•  rsj 

V  =  V 

X 

i 

< 

o 

£ 

> 

1 

1! 

•> 

o. 


Substituting  Eqn  26a  into  relations  22  through  24 
yields,  after  neglecting  second  order  terms  and  recognizing 
that  e  =  a  +  , : 


(mp  +  ma^  +  m*)  v  +  T  -  (WA  +  Wp)  =  0  (22a) 

(mp  +  ma  +  m1)v0y3  -  ma  L2(  0-+/3  )+N+(W*  +Wp)(cc+/3)  =  o  (23a) 
y  y 

(Iem  +  Ia)(«+/)  ’  L1N  =  0  d4a) 

The  normal  and  tangent  forces  (N  and  T)  are 
conventionally  represented  as: 

N  =  CN|pvp2S 

T  =  CrpIpVp  S 


where:  vp  =  +  (vy  -  u>  L2)2  ,  the  absolute  velocity  of 

the  center  of  volume  of  the  canopy. 

Following  the  assumption  of  small  oscillations, 
one  obtains: 


=  v. 


+  2  vc  v 


and  therefore: 

N  =  CN  \  p  (v02  +  2  vQ  v)TTr2  (27) 

T  =  Cip  ^  p  (v02  +  2  vQ  v)Ttr2  (28) 


where:  r  =  characteristic  radius  of  the  canopy. 

For  future  convenience,  the  related  area  of  the  coefficients 

p 

has  been  introduced  as  Tf  r  . 

Experiments  (Ref  6)  have  shown  that,  for  the 
parachutes  under  consideration,  C.p  is  approximately  constant 

over  a  relatively  large  range  of  CX  .  Substituting  relation 
28  into  Eqn  22a  yields : 

(m  +  ma  +  Tn^)v  +  CtIp(vq2  +  2v0v)Tfr2-(Wje,  +  Wp)  =  0  (22b) 

*  X 
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Reference  6  aiso  shows  that  the  normal  coefficient, 
in  the  range  of  interest,  is  proportional  to  the  angle  of 
attack.  This  may  be  expressed  as: 


ac 

(— )  a 


(29) 


where 


0  C 

( _ N\  =  slope  of  CN  versus  a  for  static  conditions 

1  aa  y  s 

Qp  =  instantaneous  angle  of  attack  of  the  canopy. 


It  should  be  noted  that  a t  Cc  .  The  angle  of 
attack  oip  of  the  canopy  is  measured  relative  to  the 

local  velocity  vector  at  the  center  of  volume  of  the  canopy. 
Thus,  when  the  canopy  oscillates  about  the  center  of  gravity 
of  the  system,  the  angle  of  attack  of  the  canopy  consists 
of  the  angle  a  (Pig  3)  and  a  contribution  aq}  induced  by 
the  rotation  of  the  canopy.  From  geometric  considerations, 
it  is  apparent  that: 


oj  L, 


CO  Lf 


a  a  = 


V 


(30) 


One  finds  oc. 


'N 


to  be: 

aP  = 

a 

+aoc 

29: 

=  ( 

0 

Cm 

0 

rv  ) 

CC  S 

cOl2 


(31) 


a+  (“  +F  )]  (32) 


Introducing  Eqns  32  and  27  into  Eqns  23a  and 
24a,  one  obtains: 


S  +  +  m*)  -  "ay  Ls  (  °  +  /?  ) 


cc 


+  (— — )  [“+  a  +/3)\iJovZitrS+(we+w  )(a+/3  )=o  (33) 

O  U  _  n  ** 


dC, 


(*«»  +  )-Ll  Q)Ja+  >] 


,  ,  +  — -  (a  +/3  )  |  |pv  2Tfr2  =  0 

0  a  sL  vo  /  ;  J  2r  o 

Equations  22b,  33,  and  34  now  represent  a  set  of 
linear  differential  equations  governing  the  motion  of  the 
parachute -load  system. 


(34) 
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iv.  DIMENSIONLESS  EQUATIONS  OP  MOTION 


It  Is  convenient  to  express  the  linearized  equa¬ 
tions  of  motion  in  a  dimensionless  form  by  introducing 
a  dimensionless  time,  distance,  mass,  moment  of  inertia 
and  velocity. 


X 


t  Vp 
r 


t  = 


L 

r 


m 

I 


m 

II  pr3 
I 

TTpr5 


(35) 


These  definitions,  substituted  in  Eqns  22b, 
and  3^+,  yield: 


+  rru )  v' 


(1  +  2v) 


+  Wp 

pvq^tt  r2 


o 


33, 


(22c) 


(mp  +  ma^  +  L2(  a"  +  /}  ") 

0C  r-  -i  w,  W 

+  L®  +  L2(  a*  +/3')\  +  - ( a  +/3 )  =  0  (33a) 

°  s  pv0  Tfr 

(lorn  +  IaK  +y3  ")  -  i  Ll  (-|^)  [«+  Lg  (  a'  +  /31)]-  0  (34a) 


where  the  prime  (')  indicates  differentiation  with  respect  to 
the  dimensionless  time  "XI 

A  further  simplification  is  introduced  by  considering 
the  definition  of  equilibrium  velocity 


Ct  I 


P 


+  W 


P 


or 


W-  +  VT 
*  P 

- 2 —  2 

Pvo  Tf  vd 


T 


(36) 
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Introducing  this  relation  into  Eqns  22c  and  33a,  one  obtains: 

(inp  +  ma^  +  rng)  v'  +  CT  v  =  0  (22d) 

(top  +  ^av  +  3  1  -  ma  L2  (  a"  +/$") 

(33b) 

+  4(-§r)s  10c+  ( a' +j}'^ + 


Rearranging  Eqns  33b  and  3^a,  one  finally  obtains 
the  new  set  of  equations: 


(mp  +  max+mi)  v»  +  CT  v  = 


0 


(37) 


m  r  1/6CN\7  rY,  rir9CN  CT 

%  l2  a  ‘ 1  v2  a  '  ^“*77v_  +  "r 


6a 


a 


(38) 


+  may  L2/3"  - 


'  .0C-N  .  - 

i  (-d^{L2  +  mP  +  ™ay+  m/ 


T 


A'  “  tA  - 


o 


(39) 


V.  FREQUENCY  EQUATION 


A.  Solution  of  the  Linearized,  Dimensionless 

Differential  Equations  of  Motion 

Upon  examination ,  it  becomes  evident  that  Eqn  37 
can  be  integrated  directly  to  give: 


CTX 


v  =  vi  e 


rap  +  may  +  ml 


(40) 


where:  v  at  X  =  0. 

One  observes,  however,  that  the  remaining  two 
equations  cannot  be  solved  as  easily  as  they  are  coupled 
together.  In  order  to  obtain  solutions  for  oc  and  fd  as 
a  function  of  X  ,  one  may  assume  solutions  of  the  form: 


9 

ii 

> 

( D 

& 

0) 

pq 

11 

(41) 

where  A,  B  and  A  are  constants.  Substitution  of  these 
relations  into  Eqns  38  and  39  yields: 

all  A  +  al2  B  = 

0 

(42) 

a21  A  +  a22  B  = 

0 

(43) 

where : 


a21  -  [xcm  +  *a]  ' 

a22  =  [Icm  +  ^a]  ^ 


ac 


N 


dcx 


) 
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Equations  42  and  43  are  homogeneous  with  respect 
to  the  constants  A  and  B.  Therefore,  the  system  will  have 
a  nontrivial  solution  only  if  the  determinant  of  the  coef¬ 
ficients  is  zero.  That  is,  if: 


D 


all  a12 
^*21  a22 


0. 


(44) 


Solving  for  D  and  setting  it  equal  to  zero,  the 
frequency  equation  of  the  system  is  obtained  as: 


(Icm  +  +  may  +  3 

2  [  -^cm  +  ^a  ^rap  +  m  t  ^  ^1  ^2  ] 


-  |  Lx (mp  +  ma 


/  3  C,T . 

+  m  )(— ^)  A 

y  a  a  s 


cT(- 


0Cn 


6a 


) 


o. 


(45) 


Equation  45  may 
a  A3 


be  written  symbolically  as: 
+  b  A2  +  c  A  +  d  =  0 


(46) 


where  a,  b,  c  and  d  are  the  coefficients  of  the  frequency 
equation  (45) . 

Before  determining  the  dynamic  stability  charac¬ 
teristics  of  the  system  by  means  of  its  frequency  equation, 
the  approach  to  the  complete  solution  shall  be  discussed. 

In  any  practical  case,  certain  initial  conditions  will  be 
known.  For  example,  at  X  =0,  one  has  a  =  a  t  y3  =  J30, 

and  (a'  +  /}')=(a ’  +/?')0*  In  general,  the  cubic 

frequency  equation  has  three  roots,  A-^,  A2  and  A^,  two  of 

which  will  usually  be  complex  numbers,  i.e.,  A  =  n  +  im. 


With  three  values  of  A  ,  one  would  have  a  general 


solution  given  by: 

AjT 

"K  2^  Agf 

(47) 

a  =  A]_e 

+  Age 

+  A^e 

Aft 

7\  2^  ^ 

(48) 

/3  =  Bje 

+  Bge 

+  B3e 
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Applying  the  previous  initial  conditions,  one  obtains: 


+  Ag  +  A3  =  oc0 

(49) 

+  B2  +  B3  =  y3o 

(50) 

(A1  +  %)  A -j_  +  (Ag  +  Bg)  A2  +  (A3  +  B3)  A3  =  (  a*  +/3')c  (51) 


So  far,  one  has  three  equations  and  six  unknowns, 
A]_,  A2,  A3,  B-j_,  Bg  and  B3.  One  must  therefore  find  three 

additional  equations  to  completely  determine  the  constants. 

These  relationships  can  be  obtained  by  utilizing 
Eqn  43.  Thus,  one  obtains  three  additional  equations  for 
the  constants  of  the  form: 


Ai 


r  dan  +  t.)  Ai‘ 


3c 


L,  L, 


(_E) 


12V  aa 


i(1om  +  A, 


/  0  CN 


-  2  B^Lg 


A' 


i  T  ^s°n 
J1 


a  t  r 

2  l(-3a  l 


(52) 


where  i  takes  on  the  values  of  1,  2,  and  3. 

One  now  has  six  equations  for  the  six  unknowns 
A^  and  B^.  The  system  is  therefore  closed  and  solvable, 

an  example  of  which  is  shown  in  Section  VI. 


B.  Stability  Criteria 

Often,  one  Is  not  confronted  with  the  problem  of 
solving  the  frequency  equation  but  one  only  wishes  to 
determine  whether  or  not  the  system  is  dynamically  stable, 
which  can  be  accomplished  by  utilizing  Routh's  criteria 
(Refs  1  and  7).  This  criteria  requires  that  for  an  os¬ 
cillating  system,  whose  oscillations  should  eventually 
approach  zero,  the  coefficient  of  the  frequency  equation 
(45)  or  (46)  must  satisfy  the  following  five  conditions: 


a  >  0 


b  >  0 


c  >  0 


d  >  0 


be  >  d 


(53) 


Examination  of  the  coefficient  of  A^  (Eqn  45) 
shows  that,  for  all  systems,  "a"  is  a  positive  term.  One 

acN 

next  observes  that  If  (-5 — )  is  negative  (i.e.,  statically 

0CX's  2 

stable  parachute)  the  coefficient  of  A  is  positive. 
Similarly,  the  coefficient  (c)  of  A  is  positive  if  the 
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Q  C 

parachute  is  statically  stable  and  d  >  0  if  ( - <  0 

v  a  a  4 

It  only  remains  to  examine  • 
means  of  substitution,  one  may  write: 

N \  /-  -  s  ~ 


0C  f 

bo-a  =  (_J.)  )  [ 


)  < 

0. 

s 

-  d. 

By 

d0N' 

)  +c 

0ot 

s 

I  +Ia+ 1 
cm  a 


The  stability  derivative  on  the  right  side  of  this  equation 
shows  again  that  a  statically  stable  parachute  is  necessary 
for  a  dynamically  stable  system.  Furthermore,  the  com¬ 
position  of  the  bracketed  term  indicates  that  the  mass  and 
inertia  terms,  the  length  dimensions  and  the  aerodynamic 

0  Q 

terms  ( - — )  and  C™  must  be  properly • balanced  in  order  to 

v9ays  1 

satisfy  this  condition  required  for  dynamic  stability. 


(54) 
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VI .  NUMERICAL  DETERMINATION  OP  THE  AMPLITUDE - 

TIME  RELATIONSHIP  OP  A  PARACHUTE  STABILIZED 
LOAD  HAVING  NEUTRAL  AERODYNAMIC  STABILITY 

✓ 

Objects  which  possess  almost  neutral  aerodynamic  * 
stability  are  frequently  decelerated  and  stabilized-  by 
means  of  an  aerodynamically  stable  parachute.  In  such 
cases,  one  desires  to  know  whether  or  not  the  system  will 
behave  with  dynamic  stability  and  how  fast  the  initial- 
oscillations  decay.  Questions  of  this  nature  can  be  answered 
through  the  solution  of  the  frequency  equation  (46) . 

For  the  purposes  of  such  a  numerical  solution, 
one  may  choose  a  ribless  guide  surface  parachute  because 
of  its  suitable  aerodynamic  stability.  For  the  determination 
of  the  apparent  mass  (ma)  one  must  know  the  enclosed  mass. 
Similarly,  the  canopy  surface  area  is  required  to  determine 
the  parachute  canopy  mass  mp.  The  Idealized  canopy  consists 
of  a  spherical  cap  and  a  truncated  cone  base.  Prom  Ref  8, 
the  volume  and  surface  area  of  the  cap  amount  to: 

Vcap  "  jV  <3R  -  hl>  (55a) 

Scap  -  2'rrHhl  •  <55b) 


Similarly,  the  volume  and  surface  area  of  the  truncated 
cone  Is : 

V  =  h0  (r2  +  rnr  +  r-j2)  (56a) 

cone  3  2  x  1  1  '  v  ' 

\ 

Scone  "  n  (r  +  rl)  [h22  +  (r  -  rl)2J  2  (56b) 

Prom  the  geometry  of  Fig  4  follows: 

o  .  JK 

h^  =  — — 2  r  r^  =  0 . 7  r 

h2  =  0.3  r  R  =  r 

Utilizing  these  values  one  obtains : 


V 


_  27  -  11  JJ^r 3 


cap 
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•V  =  0.219  Hr3 
cone 


S=ap  "  I  (3  -JT  )^rS 

(57) 

s  =  o.5i  rr  r2 

cone 
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Tho  total  volume  and  total  surface  area  of  the  canopy 
amount  to,  respectively: 

V  =  0.419  Tfr3 

S  =  1.867  Ttr2 

The  enclosed  mass  (Ref  9)  and  parachute  mass  then  become: 
me  =  pV  =  0.419  pTTr3 

m  =  crs  =  1.867  aTTr2 

AT 

where:  p  =  air  density 

a  =  mass  of  cloth  per  square  foot. 

Choosing  a  nylon  cloth  with  a  weight  of  7  oz/yd^  for  the 
canopy  material,  one  obtains  (a  =  1. 51107  x  10“^  slugs/ft  ): 


mp  =  2.821  x  10-3  TT r2 


Prom  the  definition  of  the  dimensionless  mass  m, 
finds : 


2.821  x  10~3 
pr 


one  finally 
(58) 


For  r  =  2.5  ft,  p  =  2.378  x  10  3  slugs/ft3: 


mp  =  0.475 


(59) 


Similarly,  the  dimensionless  apparent  mass  ma  can  be  derived 
from: 

ma  =  Kme  (60) 


where  K  =  0.3  from  Ref  9  for  a  nominal  cloth  porosity  of 

70  f t3/ft2-min. ,  and 

sr  =  =  0.419  (61) 

e  Tfpr3 

The  frequency  equation  (Eqn  46)  requires  the 
terms  of  the  apparent  mass  in  the  x  and  y  directions. 

As  a  first  approximation.  Ref  10  proposes  to  set: 


m 


a. 


x 
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which  amounts  to,  in  view  of  Eqns  60  and  6l: 
%x  ^  -  0-1257 


In  the  following,  one  must  know  the  center  of 
mass  of  the  system,  a  part  of  which  depends  upon  the  center 
of  mass  of  the  canopy  material.  The  calculation  is  cumber¬ 
some  but  straightforward  and  provides : 

Xm  =  0.9067  ft  (for  r  =  2.5  ft) 

mp 

Note:  Xm  is  the  location  of  the  center  of  mass 
P 

of  the  parachute  material  measured  from  the 
plane  of  the  mouth  of  the  canopy. 

Choosing  a  350  lb  load,  the  center  of  mass  of 
the  system  is  located  at  a  distance  of  0.0385  ft  above 
the  center  of  mass  of  the  suspended  weight.  The  aerodynamic 
center  of  pressure  of  the  parachute  canopy  can  be  determined 
from  conventional  three  component  measurements .  However, 
in  general  and  as  a  first  approximation,  one  may  assume  that 
the  center  of  pressure  lies  at  the  center  of  volume  of  the 
canopy.  Therefore: 

%  =  l2  =  7.34  ft 

or: 

L-^  =  L<2  =  2 . 94  . 


With  these  dimensions  and  masses,  the  moment  of  inertia  of 
the  system  is: 


=  4.36 


The  apparent  moment  of  inertia  follows  from 
Ref  11,  which  gives  Ia  for  various  canopy  shapes  as  deter¬ 
mined  by  experiment .  In  this  reference  the  apparent  moment 
of  Inertia  was  measured  about  a  point  2.66r  upstream  of 
the  plane  of  the  canopy  inlet  area.  A  dimensionless  parameter 
A  is  defined  as: 


where :  I 

a 

% 


apparent  moment  of  inertia  about  the  reference 
point 

moment  of  inertia  about  the  same  point  of  a 
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point  mass  located  at  the  center  of  volume 
of  the  canopy.  The  value  of  this  mass  is 
taken  to  be  the  mass  of  air  enclosed  in  a 
sphere  with  the  radius  of  the  inflated 
parachute . 

Thus,  It,  can  be  expressed  as: 

>,  -  § 

where:  L  =  distance  from  reference  point  to  center  of 
volume  of  the  canopy. 

Therefore,  the  apparent  moment  of  inertia  is: 

Ia  =  A  |  TTr3  p  L2  . 

The  distance  L  for  the  present  problem  is  measured  from  the 
center  of  mass  of  the  system  as: 

L  =  L2  =  2.936r. 

Thus : 

Ia  =  A  "HTp  (2. 936) 2  r3  . 

Prom  the  definition  of  Ia,  one  has: 

Ta  =  — "^15  =  t  (2.936)2  A. 

Tfpr^  5 

From  Ref  11  the  value  of  A  for  a  ribless  guide  surface 
canopy  is  0 . 187  and  one  finds : 

Ia  =  2.13  . 

Also,  since  =  350  lbs: 

m£  =  93 . 204 

dCN 

The  aerodynamic  coefficients  Cm  and  (— — )  are 

3“  s 

given  in  Ref  6  as  functions  of  effective  porosity.  For 
the  aerodynamic  coefficients  chosen  in  this  case,  the  ef¬ 
fective  porosity  amounts  to  C  =  0.025.  Thus  from  Ref  6, 
one  obtains : 


CT  =  1.08 
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6C 

(-■■w  -)  =  0.0144  per  degree 

o  a  s 

=  0.825  per  radian. 
In  this  reference,  Cm  was  defined  as: 


'm 


LN 

8 - 

3-qsr 


where  L  =  distance  from  confluence  point  to  the  apex  of 
the  canopy 

q  =  dynamic  pressure. 

With  the  definition  ~  one  finds: 


m 


_  3  L  r 

t  r 


Prom  the  geometry  of  Fig  4  follows: 


or 


CN  =  0.82  cm 


ao„  30 

-  °-82  W, 


and  finally: 


0C„ 

(-— — )  =  -O.676  per  radian. 

v  6a  ' 


The  negative  sign  has  teen  introduced  because  of  the  opposite 
sign  convention  utilized  in  Ref  12  and  the  present  report. 

Utilizing  the  previous  value  of  Crp,  one  finds: 
vQ  =  117.8  ft/sec  . 

Summarizing  all  of  the  results,  one  finds  for  a 
5  ft  diameter  ribless  guide  surface  parachute,  constructed 
of  7  oz  nylon  material  and  having  a  porosity  of  70  f t3/ft2-inin. , 
at  sea  level  conditions  with  a  350  lb  load: 


•cm 


=  4.36 
=  2.13 


(62) 
(cont . ) 
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0.419 


raa  “  ^a  -  0.1257 
ay  x 

m2  =  93.204 

Lx  =  L2  =  2.94 


-O.676  per  radian 


cT  =  1.08. 


(concl . ) 


(62) 


Substitution  of  these  values  Into  the  frequency  equation 
(45)  yields,  after  combining  terms: 

610  A3  +  276  A2  +  93.5  A  +  0.537  =  0  (63) 


Certain  terms  are  negligible  in  this  particular 
case  as  can  be  seen  by  observing  the  contribution  of  each 
term.  With  this  observation  the  original  frequency  equation 
can  be  simplified  to  the  form: 


^cm  +  ^3)  ^ 


m 1  LxL2 


ac 


m 


-r-h  (-af)7-  *  V°r 


ac 


N> 


(■ 


3a  s 

9cHx 

aa;s 


=  0 


which  provides  the  following  numerical  result: 


(64) 


605  A3  +  272  A2  +  92.6  A  +  0.537  =  0  (65) 


It  is  seen  that  the  frequency  equations  (63)  and 
(65)  are  almost  identical.  The  simplification  was  justified 
for  a  small  parachute  with  a  relatively  heavy  load.  Thus, 

If  one  were  interested  in  solving  such  a  problem,  the 
simplified  frequency  equation  (64)  could  be  used  with  very 
good  accuracy. 

Returning  now  to  the  problem  at  hand,  Eqn  63 
becomes,  on  dividing  by  610: 

A3  +  0.4525  A2  +  0.1533^  +  0.000881  =  0  .  (63a) 
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One  notices  _i'hat  all  coefficients  are  positive 
and  that  be  -  d  =  0.0684  >  0  and  thus,  this  is  a  dynamically 
stable  system. 

To  completely  specify  the  motion  of  the  system, 
the  frequency  equation  must  be  solved.  Reference  8,  page  295, 
gives  a  method  of  solving  any  cubic  equation.  Utilizing 
this  method,  one.  finds: 

=  -0.00553 

A  2  =  -0.22348  +  0.31741  i  (66) 

A3  «  -0.22348  -  0.31741  i 


Thus,  from  Eqns  47  and  48,  the  general  solutions  for  a  and 
y3  are: 


a  =  A-^e 


-0 .00553T 


/3 


Bie 


-0.00553  X 


+  AgS 


+  A3e 


+  Bge 


+  B3e 


(-0.2235  +  0.3174  i )r 
(-0.2235  -  0.3174  iyc 
(-0.2235  +  0.3174  i)tr 
(-0.2235  -  0.3174  i)r 


(67) 


(68) 


From  Eqn  52  one  obtains : 

2 


6.49  A^  +  2.92  A? 

6.49  A,2  +  2.92  A,  +  0.994 


Using  the  value  of  A ^  from  Eqn  66  yields: 


0.016307 


65.7543  -  34.7973  i 


65.7543  +  34.7973  i 


(69) 
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Using  these  relations  in  Eqn  67  one  finds: 

XfT  *2X 

Q=  0.0163076^  +  (65.7543  -  34.7973i)Bpe 

X3C  (67a) 

+  (65.7543  +  34.7973i)B3e 

As  an  initial  condition,,  outlined  in  Eqns  49 
through  51*  one  may  choose: 

at  m  =  0,  a  =  10°  =  0.1745  radians 

ft  =  0 

OC  1  +  ft'  =  0 

Thus,  one  obtains  from  Eqns  67a  and  68: 

0.01630763^  +  (65.7543  -  34.7973i)B2 

+  (65.7543  +  34.7973i)B3  =  0.1745 

+  Bg  +  B3  =  0 

and 

0.00562B1  +  (3.8749  -  28. 9648i )b2 

+  (3.8749  +  28.9648i)B3  =  0. 

Solving  these  simultaneous  equations  for  B^,  Bg,  and  B3  gives: 
Bx  =  -0.002856 

B2  =  0.001428  -  0.000191  i 

B3  =  0.001428  +  0.000191  i 

and,  from  relations  69: 

k-y  =  -0.0000466 

A2  =  0.087251  -  0.062250  i 

A3  =  0.087251  +  0.062250  i 

Using  these  relations,  one  obtains  after  a  cumbersome 
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calculation: 


a  = 


-O.OO553T 

-0.0000466e  +  0.2l434e 

-0.00553^ 

-0.002856e  +  0.002882e 


-0.22348^ 

COS  (0.31741^  -  0.62) 


-0.22348^ 

003(0.31741^-  0.133) . 


(70) 


The  values  of  CL  +  J3  are  presented  in  Table  1. 
The  amplitude -time  relationship  is  shown  in  Fig  5. 


Table  1.  CL+J3  as  a  Function  of  t 
and  t  for  the  Ribless 
Guide  Surface  Parachute 


T 

t 

a  *fi 

(sec . ) 

(deg. ) 

0 

0 

10.000 

2 

0.042 

7.780 

4 

0.085 

2.914 

6 

0.127 

0.736 

8 

0.170 

-0.884 

10 

0.212 

-1 . 384 

12 

0.254 

-1.002 

14 

0.297 

-0 . 569 

16 

0.339 

-0.237 

18 

0.382 

-0.066 

20 

0.424 

-0.036 

22 

0.466 

-0.056 

24 

0.509 

-0.102 

26 

0.551 

-0.136 

28 

0.594 

-0.152 

_ 20 _ 

0.636 

-0.148 
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Fig  5,  QC+p  as  a  Function  of  Time  for  the  Ribless  Guide 
Surface  Parachute 
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